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flj ' Abstract 

T^ We study an extension of the Higgs sector of the Minimal Supersymmetric Standard 

Model (MSSM), considering the effects of new degrees of freedom at the TeV scale, 
^NJ and allowing for sources of CP violation beyond the MSSM (BMSSM). We analyze the 

J> impact of the BMSSM sources of CP violation on the Higgs collider phenomenology 

'^ and on low energy flavor and CP violating observables. We identify distinct Higgs 

__. collider signatures that cannot be realized, either in the case without CP violating 

fK^ phases or in the CP violating MSSM, and investigate the prospects to probe them at 

the Tevatron and the LHC. The most striking benchmark scenario has three neutral 
Higgs bosons that all decay dominantly into W boson pairs and that are well within 
the reach of the 7 TeV LHC run. On the other hand, we also present scenarios with 
three Higgs bosons that have masses M//. > 150 GeV and decay dominantly into 
bb. Such scenarios are much more challenging to probe and can even lie completely 
outside the reach of the 7 TeV LHC run. We explore complementary scenarios with 
/\ standard MSSM Higgs signals that allow to accommodate a non-standard Bg mixing 

phase as indicated by DO, as well as the excess in Bg -^ M^/^^ candidates recently 
reported by CDF. We find that, in contrast to the MSSM, a minimal flavor violating 
soft sector is sufficient to generate significant corrections to CP violating observables 
in meson mixing, compatible with EDM constraints. In particular, a Bg mixing phase 
Sitxp ^ 0.15, can be achieved for specific regions of parameter space, compatible with 
all the presently available experimental constraints on fiavor observables. Such a non- 
standard Bs mixing phase would unambiguously imply a sizable suppression of S^Ks 
with respect to the SM prediction and a BK{Bs -^ fi^fi~) close to its 95% C.L. upper 
bound of 1.1 X 10~*^. 
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1 Introduction 

Extensions of the Standard Model (SM) that are based on Supersymmetry (SUSY) are 
arguably the most popular models of New Physics (NP) at the TeV scale. In particular 
the Minimal Supersymmetric Standard Model (MSSM) is one of the most studied NP 
models. On the other hand, arguments based on the naturalness criterion suggest that 
the minimal Higgs sector of the MSSM should be extended in order to raise the tree level 
mass of the lightest Higgs boson above its MSSM limit of M^ cos 2/3 [1, 2, 3]. Many such 



Supersymmetric models that contain physics Beyond the MSSM (e.g. an enlarged gauge 
sector, additional singlets) have been constructed (see e.g. [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 
14]). As long as the scale of the beyond MSSM physics M is sufficiently larger than the 
mass scales of the MSSM Higgs sector, it is possible to describe the effects of the new 
degrees of freedom in an effective theory approach, where higher dimensional operators are 
added to the MSSM [15, 16, 17, 18, 19, 20, 21]. It was shown that in such an effective field 
theory approach (the so called BMSSM), the lightest Higgs boson can easily be enhanced 
at the tree level and reach a mass of around 200 GeV, as long as M is not larger than a 
few TeV. 

In this work we consider the BMSSM framework where the MSSM Higgs sector is ex- 
tended by the leading SUSY preserving and SUSY breaking dimension 5 operators [16]. Ex- 
tensions of the MSSM including also dimension 6 operators have been studied in [19, 20, 21] 
and examples of possible UV completions have been presented in [19]. One interesting fea- 
ture of the model in [16] is, that in contrast to the MSSM it allows for CP violation in 
the Higgs sector already at the tree level. Most studies of this framework in the literature 
assume the absence of new sources of CP violation and consist of analyses of the vacuum 
structure of the model [22, 23], dark matter [24, 25, 26] and the Higgs collider phenomenol- 
ogy [27]. The possible effects of CP violation induced by the higher dimensional operators 
have been mainly studied in the context of electro- weak baryogenesis [28, 24, 29]. 

Differently, in the first part of this work we study the impact of the BMSSM sources 
of CP violation on the Higgs collider phenomenology, extending the analysis of the CP 
conserving case of [27]. In [27] the expected signals of the CP conserving BMSSM at the 
Tevatron and at the LHC have been worked out in detail. Higgs production and decay 
patterns that are markedly different from the MSSM have been identified. One example are 
scenarios where both CP-even scalar Higgs bosons decay dominantly into a pair of gauge 
bosons. In this work we identify characteristic collider signatures of the BMSSM with CP 
violation that cannot be realized, either in the case without CP violating phases or in the 
CP violating MSSM and investigate the prospects to probe the model at the Tevatron and 
the LHC. 

In the second part of this work we analyze possible characteristic signals of the higher 
dimensional operators in flavor physics. As the higher dimensional operators mainly modify 
the spectrum and couplings of the neutral Higgs bosons of the MSSM, significant deviations 
from the MSSM predictions are expected in those flavor observables that are highly sensitive 
to the exchange of neutral Higgs bosons. Consequently, we analyze the rare B^is — )■ fJ.^fJ-~ 
decays as well as B^^s — B^^s mixing, that can receive sizable NP contributions in the large 
tan/3 regime from Higgs and double Higgs penguins, respectively. In particular, we study 
if the BMSSM with Minimal Flavor Violation (MFV), i.e. the BMSSM with no additional 
sources of flavor violation apart from the CKM matrix [30, 31, 32] but new sources of CP 
violation from the dimension 5 operators, can accommodate a large phase in Bg mixing 
as indicated by recent experimental results from Tevatron, especially the very recent DO 
result on the like-sign dimuon charge asymmetry [33]. In view of the excess of B^ — )• /U^/U^ 



candidates recently reported by CDF [34], we in particular also analyze the correlation 
between the Bs mixing phase and the BR(i?s -^ A*^^") in the BMSSM with MFV. A 
complementary discussion of other flavor observables in the BMSSM appeared recently 
in [35]. 

The paper is organized as follows. In Sec. 2 we review the extended Higgs sector of 
the BMSSM with sources of CP violation at tree level. We discuss the Higgs potential in 
presence of the higher dimensional operators and the vacuum structure of the model. We 
analyze the Higgs spectrum emphasizing the possible role of the new CP violating phases. 
In Sec. 3 we discuss constraints coming from Electric Dipole Moments (EDMs) that are in- 
duced by the new phases appearing in the Higgs sector. The Higgs collider phenomenology 
of the model is discussed in section 4. We address constraints coming from direct Higgs 
searches at LEP and Tevatron and present benchmark scenarios showing Higgs spectra 
and couplings that are specific to the BMSSM with CP violation. We outline the most 
promising Higgs search strategies in these scenarios. In Sec. 5 we discuss distinct BMSSM 
signals in flavor physics, concentrating on the phase of Bg mixing and its correlation with 
the rare Bg — ;■ /^^/U^ decay. We conclude in Sec. 6. The appendices contain some details 
on the chargino, neutralino and squark masses in the BMSSM, a short discussion about 
electroweak precision tests, as well as a compendium of loop functions. 

2 The Extended Higgs Sector of the Model 

We study the framework first presented in [16], where the leading higher dimensional 
operators are added to the MSSM Higgs sector. The scale M at which these operators 
arise is assumed to be not far above the TeV scale. 

2.1 The Higgs Potential 

At the 1/M order, the most general Higgs superpotential reads [16] 

W = fiHuHd + ^[HuHdy , (1) 

where H^ and Hd are the Higgs superfields with hypercharge +1/2 and —1/2, respectively 
and we denote HuHd = H^H^ — H^H^. The dimensionless parameter uj is taken to be of 
order 1 and complex. 

In addition to the 1/M suppressed term in the superpotential, a corresponding SUSY 
breaking term is added to the Lagrangian [16], 

CDa'^iH^Hdf, (2) 

where a is another free parameter of order 1 and complex. The scale ms is the scale of 
the SUSY breaking terms of the physics beyond the MSSM and therefore ms <C M is 



necessary to integrate out the complete SUSY multiplets of the new degrees of freedom at 
a common scale M. The scale ras can be connected to the scale of the soft SUSY breaking 
terms of the MSSM (i.e. the sfermion and gaugino masses), but it is possible to allow for 
sizable differences between these scales. To ensure validity of the effective description of 
the BMSSM physics in terms of higher dimensional operators, mg (as well as jjl) has to be 
sufficiently small compared to Af , for a and uj of 0(1). 

At the renormalizable level, the resulting tree level scalar potential then reads 

= (^k + \f^\^)HlHu + {ml^ + \^l\^)H\Hd + {Bii{H^Hd) + h.c.) 



+ [^^iH^H,f-^{H^H,){HlH^ + HlH,) + h.cj , (3) 



and for later convenience we define 

A, = |W*=.^, A« = |A**^i^. (4) 

The 1/M operator in the superpotential leads to two additional non- renormalizable dimen- 
sion six terms 

Vg = ^^{HuHd){HiH\){HiHu) + ^{HuHa){HiHl){H\Hd) , (5) 

with As = Ag = |wp. These terms are essential to stabilize the Supersymmetric electro- 
weak Symmetry Breaking (sEWSB) vacuae analyzed in [22]. Thanks to these non-renor- 
malizable terms, the potential is automatically bounded from below. 

At the I/M"^ order, there can be additional operators in the Kahler potential that mod- 
ify the quartic couplings of the Higgs potential. Their possible impact has been analyzed 
in [19, 20]. In this work instead, we focus on the leading effects generated by the 1/M op- 
erators in (1) and (2). While additional 1/M^ operators can lead to an additional increase 
of the lightest Higgs mass, we do not expect them to change our main conclusions. 

The three parameters B^, a and w can in general be complex. We follow the usual 
convention adopted in studies of the MSSM and absorb the phase of B^ by a rephasing 
of the two Higgs doublets. In addition, we will assume all possible complex parameters of 
the MSSM (e.g. gaugino masses, /i parameter etc.) to be real. If these parameters were 
complex, CP violating effects in the Higgs sector would be possible at the 1-loop level and 
would add to the tree level effects. In this work we consider the phases of a and w as the 
only Beyond the Standard Model sources of CP violation. 



2.2 The Minimum of the Potential 

We parametrize the Higgs fields as^ 

Hu = e^'- f w f f ^ . ^ . H, = e^^^ f A^"'^ + ^^ + ''''^\ , (6) 

\^[Vu + K + iau)l ' " \ R- J ' ^ ^ 

where Vu = v sin /3 = us^ and v^ = v cos /3 = i;c^ with t; = 246 GeV are the two vacuum 
expectation values (VEVs) and tan/3 = tp = Vu/vd is their ratio. 

While the relative phase of the Higgs fields 9u — 0d can be rotated away by a U{1)y 
transformation, 6 = 9^ + 9^, is a physical phase. Therefore, the following three extremal 
point conditions have to be satisfied at the minimum of the potential 

dV dV dV ^ ^^^ 



dReHu dReHd 89 

Using the first two conditions in Eq. (7), the two soft masses mn^ and mH^ can be traded 
for V and tan/3. The third condition determines the phase of the Higgs VEV, 0, as a 
function of the phases of a and uo and reads 

v^C/jS/jIAsI sin((/)5 + 29) + v^\\fs\ sin((?!)6 + 9)- 2Bfism9 = . (8) 

Contrary to the MSSM, the BMSSM predicts in general already at tree level a non-zero 
phase of the Higgs fields at the minimum. In Sec. 3, we will show that constraints coming 
from Electric Dipole Moments (EDMs) imply that the phase of the VEV is typically rather 
small 9 < O(O.l). 

To ensure that the stationary point of the potential defined by (7) is a minimum, 
the corresponding Hessian has to be positive definite. However, the conditions to have a 
stationary point and the requirement on the Hessian do not necessarily lead to a unique 
solution. If the quartic couplings along the D-flat direction are negative, then a second 
minimum in the Vu — Vd plane may arise for large field values that is stabilized by the 
1/M^ suppressed terms in (5). If this second minimum is deeper, the physical minimum 
at V = 246 GeV is unstable and will decay into the second deeper minimum. In Fig. 1 we 
show in the ^ — ms plane the region that is allowed by the requirement of absolute vacuum 
stability. The criterion of absolute stability particularly leads to stringent upper bounds on 
the value of the /i parameter. In accordance with the findings of [23], we observe that the 
portion of allowed parameter space shrinks for smaller values of the charged Higgs mass. 
This behavior with Mfj± holds, varying the value of the NP scale M or the two phases 
of CO and a, even if their value can have a rather sizable impact on the allowed values for 
yU. Combining the requirement of vacuum stability with the bound on the mass of the 
chargino, that leads to a lower bound on the absolute value of ^ of ~ 100 GeV, the allowed 



^Here we neglect the possibility of charge breaking vevs. As shown in [23] this is a good approximation 
in regions of parameter space that lead to a stable vacuum. 
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Figure 1: Constraints in the /i — ms plane from vacuum stability. The region within the blue, 
(green, red) contour leads to an absolute stable electroweak vacuum for a charged Higgs mass of 
150 GeV (250 GeV, 350 GeV). The NP scale M is fixed to 2 TeV, tan/3 = 2, |w| = |a| = 1, 
m^ = 800 GeV and At — 2m(. In the left plot both a and uj are real and positive, while the right 
plot shows the situation with maximal phases Arg(Q;) — Arg(w) — tt/2. The gray band is excluded 
by direct bounds on the chargino mass. 

regions for fi are strongly constrained, especially for low values of the charged Higgs mass 
(see the blue contours in Fig. 1 corresponding to Mjj± = 150 GeV). 

However, one should keep in mind that the requirement of absolute vacuum stability is 
rather conservative: it would be sufficient to impose that the EW vacuum is metastable, 
provided that its life time is longer than the age of the universe. This possibility has 
been discussed in [23] in the context of the BMSSM without CP violation. Using simple 
analytic approximations for the bounce action [36], we checked that the viable parameter 
space indeed opens up to some extent if we allow for a metastable vacuum. However, to 
be conservative, we require always absolute stability of the vacuum in the discussion of the 
Higgs phenomenology. 



2.3 The Higgs Spectrum 

We now briefly review the effects of the higher dimensional operators introduced in the 
previous section and in particular of the new physics phases on the Higgs spectrum. In 
order to keep a clear comparison to the case without CP violation, we write the mass 
matrices of the neutral Higgs bosons in the basis that would diagonalize it in absence of 
new sources of CP violation 



S/3 -Ci3 

C/3 8/3 J \ad 



(9) 



The angle a is given by 




sin 2a = 


Ml + Ml . ^^ , 2t;2|A6|cos((/.6 + 20) 
Ml-Ur'^^^' Ml -Ml 


cos 2a = 


Ml -Ml v''\\^\cos{<l^^ + 2e) 
Ml-Mf""-^ Ml- Ml ' 



(10) 
COS 2/3. (11) 

In absence of CP violation, Ma is the mass of the pseudoscalar Higgs^ 

MaS0Ci3 = Bfi COS 6 — — jAel cos(06 + 6*) — f |A5|s/3C/3Cos((/>5 + 26*) , (12) 

and Ml and Ml are the masses of the two scalars, namely the eigenvalues of the mass 
matrix 



c 



2 „ „ \ / „2 



^^""^^ ^ m2 f ^P -C/35/3A 



Ml = Ml[ ^/ :r +M|. ^ , 

+t;2|A6|cos(06 + ^)PT'' o.LV^'l^5|cos((/>5 + 2^)ff °)-(13) 



y 1 2cfjSi3j y s 



While the Goldstone boson G is not affected by the presence of CP violation, the 3 physical 
Higgs bosons mix in presence of CP violating phases. In the basis {h,H,A), their mass 
matrix can be written as 

/Ml Ml A 
Mj,= { Ml Ml A , (14) 

\MIa Ml^ Ml J 

where the mixing terms are given by 

MIa = -y(C;3+«|A5|sin((/)5 + 20)-2S/3_„|A6|sin((/)6+^)) , (15) 

^2 
MIa = -y (s^+a|A5|sin((/.5 + 2^)-2C/3_„|A6|sin(06+^)) . (16) 

The Higgs mass matrix (14) can be diagonalized by an orthogonal matrix O 

O^MlO = di8.g{Ml^,Ml^,Ml^) , (17) 

where Ml_ are the three eigenvalues. 



^ In the CP violating scenario we are studying, Ma as well as Mh and Mh are only auxiliary parameters 
and not physical masses. 
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Figure 2: The lightest neutral Higgs mass in the Arg(a) - Arg(a;) plane for two different values 
of tan/3 — 2, 5. The remaining parameters are fixed to: |a| = |w| = 1, /i = ms = 150 GeV, 
M = 1.5 TeV, Mn± = 200 GeV, mi = 800 GeV, At = 2mi. In the red hatched region the 
electro-weak symmetry breaking vacuum is only a local minimum of the potential. 



In order to get an analytical understanding of the dependence of the Higgs spectrum 
of the model on the several phases, we give approximate expressions for the masses in the 
decoupling limit Ma ^ Mz, performing a simultaneous expansion in l/tan/3 and 1/M 
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(18) 



M 



~ MI + -IA5I cos((/<5 + 2^)-l , 



H2 



Mfj^ ~ Mi + -|A5|(cos(05 + 2^) + l) . 



(19) 
(20) 



In the second line of (18) we also included the dominant 1-loop corrections to the lightest 
Higgs mass for which they are most relevant. As one notes from the approximate tree 
level part of the expression in (18), the BMSSM effects on the lightest Higgs mass are 
relevant for not too large values of tan /3 and of the NP scale M (entering the approximate 
expression through Ag). Their sign depends mainly on the phase of \q, namely the phase 
of u. Dependence on the phase of A5 and correspondingly on the phase of au arises first at 
the order l/(tan^/3M), i.e. it is only relevant for very small values of tan/3. This feature 
is also illustrated in Fig. 2 where we compare the value of the mass of the lightest Higgs 



boson in the Arg(a) — Arg(ct;) plane, for tan/3 = 2 (left) and tan/3 = 5 (right), fixing the 
remaining free parameters of the model to a reference point. 

In Fig. 2, as well as in the remainder of this work, in the numerical computation of 
the Higgs spectrum, we always include the 2-loop supersyminetric corrections to the Higgs 
potential as given in [37, 38] and diagonalize the 3x3 Higgs mass matrix (14) numerically. 

The two plots in Fig. 2 show the decoupling of the NP effects on the lightest Higgs mass 
with tan/3. For tan/3 > 10, the lightest Higgs mass differs from the MSSM expectation 
only by few GeV. From the figure, it is also evident that the maximal values for the lightest 
Higgs mass are obtained in the CP conserving case Arg(a) = Arg(a;) = 0. 

For completeness, we also give the analytical expression for the charged Higgs mass 

Ml± ~ Ml + M^ + y|A5|cos(,^5 + 20) 

~ BixcosO ^lAel cos((/)6 + 0) - ^jAsI cos{(t)^ + 29) + M^ . (21) 

cpsp cpsp 1 1 

In presence of CP violation it is customary to characterize the Higgs sector in terms of 
tan/3 and Mu± instead of tan/3 and Ma, given that Ma is no longer a physical mass. 
Unlike in the MSSM however, we note that in the BMSSM there is not necessarily a one 
to one correspondence between Bjjl and Mjj±, already in the CP conserving case. In the 
upper plot of Fig. 3 we show in an example of a CP conserving scenario the charged Higgs 
mass as a function of Bfi for several values of tan/3. Choosing for example tan/3 = 20, 
we observe that charged Higgs masses between 100 GeV and 250 GeV can be realized by 
two different choices of Bfj,. The non-monotonic dependence of Mjj± on Bfi arises because 
the phase of the Higgs VEV that enters Eq. (21) changes by varying Bfj,. In fact, in the 
CP conserving case, for large enough values of Bfi, the minimization condition (8) implies 
6 = and the charged Higgs mass decreases with decreasing Bfj, (see solid curves in Fig. 3). 
On the other hand, for very small values of B^ and a positive Ag, the condition (8) implies 
6 = TT and the charged Higgs mass increases again for further decreasing B^ (see dashed 
curves in Fig. 3).^ The range of charged Higgs masses that can be realized by two different 
values of B/j, becomes smaller for smaller tan/3 and \Xq\ and eventually vanishes. 

Given the fact that one Higgs mass can potentially be realized by two different values 
of -B/i, we conclude that fixing the charged Higgs mass, tan/3 as well as A5 and Ag does 
not uniquely determine the Higgs sector of the theory. This is further illustrated in the 
lower plots of Fig. 3 that show the neutral Higgs spectrum as function of the charged Higgs 
mass with all other parameters fixed. The solid curves correspond to large values of the 
Lagrangian parameter Bfi (and consequently 9 = 0), while the dashed curves correspond to 
small B^ (and consequently 9 = tt). In particular for the lightest Higgs mass one observes 
a significant shift between the two cases. 



■^For intermediate Bfi there exists a region of spontaneous CP violation where the phase of the Higgs 
VEV changes continuously from to tt. That region is however not phenomenologically viable as the neutral 
Higgs spectrum becomes tachionic. 
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Figure 3: Top: The charged Higgs mass Mh± as a function of the B/i parameter for different 
values of tan /3 as indicated in the plot. Bottom: the three neutral Higgs masses M//. as a function 
of Mii± for tan/3 = 20. In all plots the remaining parameters are fixed as: \o\ = \uj\ = 1, 
^ = TO5 = 150 GeV, M = 1.5 TeV, m^ = 800 GeV, At = 2mi, Arg(a) =Arg(w) = 0. 

In the scenarios that we are analyzing in the remainder of this work, tan f3 is fixed to 
a rather small value and we checked that the given charged Higgs masses do fix the Higgs 
sector of the model and specifying the Bfj, term is not necessary. 
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Figure 4: Contributions to the electric and chromoelectric dipole moments most relevant for 
our analysis. Diagram a) is a 1-loop Wino-Higgsino-sneutrino contribution to the electron EDM, 
while diagram b) is a gluino-down squark contribution to the down quark (C)EDM. Diagram c) 
represents the dominant 2-loop Barr-Zee type contribution to both the electron EDM and down 
quark (C)EDM. 

3 Constraints from Electric Dipole Moments 

In this section we discuss the constraints coming from Electric Dipole Moments (EDMs) 
on the new CP violating phases arising in the Higgs sector. 

EDMs are known to be highly sensitive probes of CP violation in extensions of the 
Standard Model [39]. Indeed the current bounds on experimentally accessible EDMs like 
the ones of Thallium, Mercury and the neutron [40, 41, 42], 



dxi < 9.4 X 10^2^ ecm @ 90% C.L. , 

dug < 3.1xl0~2%cm @ 95% C.L. , 

dn < 2.9 X lO"^*^ e 



cm 



C.L. , 



(22) 
(23) 
(24) 



generically lead to very tight constraints on new sources of CP violation that can be present 
in extensions of the Standard Model. In particular in the MSSM with SUSY particles at 
the TeV scale, the flavor diagonal CP violating phases of e.g. the gaugino masses, the 
higgsino mass parameter and the trilinear couplings are strongly constrained [43, 44, 45]. 
In the following, we take all MSSM parameters real and concentrate on effects of the 
phases of the BMSSM parameters a and oj. In this setup important NP effects on the 
experimentally accessible EDMs are induced by contributions to the lepton and quark 
(chromo)electric dipole moments ((C)EDMs) both at the 1-loop and at the 2-loop level. ^ 
Concerning the 1-loop contributions, note that the higher dimensional operators do not 
only modify the Higgs sector of the model, but, after electro- weak symmetry breaking, 
also lead to complex entries in the chargino, neutralino and squark mass matrices (see 
appendix A). In particular, the most important effect arises from the phase of the Higgs 



''We explicitly checked that in the scenarios that we consider in the present work, four fermion opera- 
tors [46] are always sub-dominant. Still, they are consistently included in our numerical analysis. 
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VEV that enters these mass matrices, resulting in potentiaUy large 1-loop contributions 
to both the electron and quark (C)EDMs even if the gaugino and Higgsino masses as well 
as the trilinear couplings are assumed to be real. In our numerical analysis we include the 
full set of SUSY 1-loop contributions calculated in the mass eigenstate basis following [43] . 
The most important diagrams are a Higgsino- Wino-sneutrino loop for the electron EDM 
and a gluino-down squark loop for the down quark (C)EDM. They are shown in Fig. 4 and 
have the following approximate expressions 
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(26) 



These expressions show clearly that the 1-loop EDMs are mainly induced by the Higgs 
phase 6. The loop functions fd, fd and fe depend on the mass ratios Xg = M?/rh'^, 
x^ = ij? /m? and X2 = M2/m? and their analytical expressions are given in appendix C. 
For simplicity we set all sfermion masses to a common mass scale m in the approximate 
expressions we show in this work. The e terms appearing in the all order tan j3 resummation 
factors, arise from non-holomorphic corrections to the down quark and electron Yukawa 
couplings. Such corrections become relevant in the large tan/3 regime and read 
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(27) 



Here we only included gluino, Higgsino and Wino loops, but neglected Bino loops that are 
typically not relevant. The loop functions /i and /2 can again be found in appendix C. 

It is interesting to note that the 1-loop contributions that are sensitive to the phase 
of the Higgs VEV are those that involve non-holomorphic couplings of the electron and 
down quark to the up-type Higgs and are therefore tan/3 enhanced.^ Similarly, also in the 
expressions of the tan/3 resummation factors the phase of the Higgs VEV appears. 

At the 2-loop level the most important contributions come from Barr-Zee diagrams, includ- 
ing a top quark loop, that are directly sensitive to the mixing of the scalar and pseudoscalar 
Higgs bosons [47]. The diagrams are again shown in Fig. 4. Their dominant contribution 



^One loop contributions to the up quark (C)EDM that are sensitive to the phase of the Higgs VEV are 
thus l/tan/3 suppressed and only relevant for a very small tan/? ~ 1. 
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is tan/3 enhanced and can be approximated by the following expressions 



''""/' = i^r-"" 



tan/3 






2 



X 



E^03,f^02. + ^Oi.)/(-;^) , (28) 

\ '^ ' ' <i } me Re[l + eitf}] " ^ \ 3 ' ae„^ 8 j ^ ^ 

Here, 0,^ are the elements of the matrix O defined in (17) that diagonalizes the neutral 
Higgs mass matrix. Subleading contributions that are not enhanced by tan/3 can become 
important for small values of tan/3. The 2-loop function / that enters Eq. (28) can be 
found in appendix C. The expression (28) can be further expanded in the decoupling limit 
and performing an expansion in 1/M. We find 

X ^^ I ^{y ^^ Cry ^-, \ n I TTi-t 






^^h 



^^|A5|sin((/>5 + 2g) v^\\^l^{2^^ + 2e)\ 

with the mass ratio z = m1/M\. Eq. (30) clearly shows that the 2-loop EDMs are directly 
induced by the new phases of the higher dimensional operators. In the scenarios we con- 
sider, the Higgs bosons are always light compared to sfermions and the 2-loop contributions 
can compete with or even dominate the 1-loop contributions discussed above. If the stop 
mass is small, then in addition also 2-loop diagrams with stop loops are often relevant. For 
large values of tan /3, (s)bottom and (s)tau loops can also become important. Even though 
the sparticle masses are rather large in the scenarios that we consider in the following, 
we include the full set of 2-loop Barr-Zee contributions from [48, 47, 46] in our numerical 
analysis. 

Expressing the experimentally accessible EDMs of Thallium, Mercury and the neutron 
through the quark and electron (C)EDMs induces sizable uncertainties related to QCD, 
nuclear and atomic interactions. Approximately one finds the following relations [39, 49] 

dxi ^ -5854 , (31) 

^Hg ^ 7 X m-^e{du - dd) + IQ-^de , (32) 

dn ^ 1.4((id-0.25(i„) + l.le((id + 0.5(i„) . (33) 

The quark (C)EDMs in the above expressions are understood to be evaluated at a scale of 
1 GeV. Expressions for the running of the EDMs down from the high matching scale can 
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Figure 5: Bounds from the EDMs in the Arg(Q:) - Arg(aj) plane for two different values of \a\ — 
1, 0.5. The remaining parameters are fixed to tan/3 — 2, |a;| = 1, fi — ms — 150 GeV, M = 1.5 TeV, 
Mh± = 200 GeV, m = 800 GeV, At = 2m, Ab = A^ = and Mg = ilvh = 6Mi = 1200 GeV. 
The solid and dashed black lines correspond to the Thallium and Mercury EDMs respectively. 
The neutron EDM does not lead to constraints in the considered scenarios. The dotted blue lines 
indicate the values of the phase of the Higgs VEV 9. In the red hatched region the electro-weak 
symmetry breaking vacuum is only a local minimum of the Higgs potential. 



be found e.g. in [50]. While the prediction for the Thallium EDM is rather robust, the 
uncertainty in the neutron EDM is estimated to be at the level of 50% and the expression 
for the Mercury EDM is only accurate up to a factor of 2-3 [49]. We take these uncertainties 
into account when evaluating the corresponding constraints. 

As the dominant contributions to the EDMs are tan (3 enhanced, larger values of tan /3 
lead to stronger constraints. Here we restrict the discussion to the low tan/3 regime, where 
we expect the most interesting Higgs collider phenomenology (see Sec. 4). A detailed 
treatment of EDMs for large tan/3 will be important when we explore a complementary 
region of parameter space analyzing interesting effects in the flavor phenomenology (see 
Sec. 5). 

In Fig. 5 we show the constraints coming from the EDMs in the Arg(6t;) — Arg(a) plane 
for two example scenarios with tan/3 = 2. We observe that the most stringent constraint 
comes from the Thallium EDM that is dominantly induced by the 1-loop chargino con- 
tribution to the electron EDM. From (26) one expects that the allowed region therefore 
corresponds to small values of the phase of the Higgs VEV as it is also shown in Fig. 5, 
where the values of 9 are indicated by the dotted blue contours. In (26) we neglected 
additional 1/M suppressed corrections that can be incorporated by replacing the Higgsino 
mass with an effective term /ie*^ — t- /ie*^ — ojj^spcpe^'^^ . Therefore, the region compatible 
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with the ThalHum EDM is shghtly tilted with respect to the ^ = Hne in the plots of 
Fig. 5. 

Concerning the Mercury EDM, we find that generically it is mainly induced by the 
2-loop CEDM of the down quark in Eq. (29). Interestingly, also the regions allowed by the 
Mercury EDM correspond to small values of the phase of the Higgs VEV in Fig. 5. Given 
its large theory uncertainty however, the Mercury EDM is typically less constraining than 
the Thalhum EDM. 

It is instructive to derive from Eq. (8) an analytical expression for the 9 = line 

where the EDM constraints are in general easier satisfied. Defining cp^ = Arg(a;) and 

(pa = Arg(a), one has 

• A {alms . , . , . 

sm(pai = -spcp sm((/)Q + 0^) . (34) 

^^ 

This shows that for smaller values of |a| and larger values of tan/3, the allowed region 
corresponds to smaller values of (/>^. Usually one finds that the full range — tt < Arg(a) < vr 
remains allowed if Arg((^) is chosen appropriately. For this reason, in the following we will 
study observables in the Higgs collider phenomenology mainly as function of Arg(a). 

4 Higgs Collider Phenomenology 

4.1 LEP and Tevatron Bounds 

In this section we study the LEP and Tevatron bounds on the Higgs bosons of the model. 
In order to evaluate the bounds from direct Higgs searches on the BMSSM parameter 
space, knowledge of the couplings of the Higgs bosons is required. It is convenient to work 
with effective couplings that are normalized to the corresponding SM couplings. For the 
effective HiZZ and HiWW couplings, that are essential both for the Higgs production at 
LEP and the Higgs decays at the Tevatron and the LHC, one finds expressions that have 
the same structure as in the MSSM with CP violation 

izZHi = iwWHi = Sp-aOli + Cp_a02i ■ (35) 

We stress however that the angle a and the rotation matrix O differ from the MSSM. 
They are highly sensitive to the higher dimensional operators of the BMSSM and given in 
Eqs. (10), (11) and (17). The couplings obey the relation 



T.' 



izzH. = 1 ■ (36) 

Possible deviations from the above expressions arise from 1/M^ suppressed operators [19] 
that we however do not consider here. Similarly, in the considered framework also the 
Higgs-quark and Higgs-lepton couplings have the same structure as in the MSSM with 
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CP violation, with the angle a and the rotation matrix O given by their BMSSM expres- 
sions (10), (11) and (17). 

We calculate the effective Higg and Hi^^ couplings from the ratios of the LO decay 
widths to the LO SM decay widths 

^^^''^ r(/7,^77)^S' ^ ' 

2 ^ rjH, ^ gg)^0 ^ ajgg ^ Hi) 

^99H. r(/7, ^ gg)\o^ a{gg ^ //.)sM ' ^ ^ 

In our calculation we incorporate the full set of SM and SUSY particle contributions using 
the expressions in [51]. As stated in Eq. (38), we assume that the effective H^gg coupling 
approximates the ratio between the gg — t- Hi production cross sections at Tevatron and 
LHC and the corresponding SM cross sections. This approach has also been adopted in [27], 
where it has been explicitly checked that it leads to results in the CP conserving BMSSM 
that are accurate within 5% - 20%, depending on tan/3. We do not expect this to change 
in presence of CP violation. 

To check compatibility with Higgs searches at LEP and Tevatron, we use the latest 
version of Higgsbounds [52, 53] in the effective coupling approximation. To obtain the 
total decay width of the Higgs bosons of the BMSSM, we make use of the results collected 
in [51], replacing the MSSM couplings with the BMSSM ones. In the computation of the 
several partial decay width, we use ratios of partial decay widths in the BMSSM and the 
SM and multiply the results with the state of the art SM partial decay widths obtained 
from HDECAY [54], when applicable. We also cross checked our results using a version of 
CPsuperH [51, 55] with appropriately modified couplings. 

We remark that by default, Higgsbounds uses the latest combined SM Higgs exclusion 
from the Tevatron [56] only for Higgs bosons that satisfy very restrictive requirements on 
their "SM-likeness" . In particular, a Higgs boson is considered SM-like by Higgsbounds 
if its different cross sections normalized to the SM values differ at most by 2% from a 
common scale factor. For Higgs bosons that do not satisfy this requirement, the strongest 
constraint used is then typically the gg ^ Hi ^ WW analysis in [57] that is considerably 
weaker. In our BMSSM scenarios discussed below, the most distinct cases are those in 
which the three neutral Higgs bosons share couplings to gauge bosons and fermions in a 
non SM-like way. In order to get a reasonable estimate of the current Tevatron bounds on 
the BMSSM parameter space, we therefore consider in addition to Higgsbounds also the 
latest Tevatron exclusion. As in the region that is excluded by Tevatron Higgs searches, 
the dominant process is gg ^ Hi ^ WW, we compute the corresponding cross section in 
the BMSSM normalized to the SM and apply the bounds given in [56] for Higgs masses 
in the range 145 GeV < Af//. < 200 GeV, where the effect of vector boson fusion and 
associated production is minimal. 

Also in the low mass region, Higgsbounds does not use the combined Tevatron exclu- 
sion [58] , but applies separately the different analyses entering the combination, for most of 
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Figure 6: Lightest Higgs boson mass as a function of the charged Higgs boson mass (left) and the 
phase of a (right). Shown are the points excluded by EDM constraints (in blue), LEP (in orange) 
and Tevatron (in red). Green points are allowed. The scan is performed fixing tan/3 = 2, mg = 
^ = 150 GeV, a common squark mass of m — 800 GeV, a common slepton mass of rhi = 1100 GeV, 
trilinear couplings At = 2m, Ai, = At = and gaugino masses Ah = 3M2 = 6M1 = 1200 GeV 
and scanning over a, oj, M, Mh± in the ranges a — (0.5 — 2)e*^'^s(") with Arg(a) g [0, 2n], oj = 
(0.5 - 2)e-5^''s("), M = (1 - 3) TeV and Mh± < 350 GeV, respectively. 

the BMSSM parameter space. We find however that even the combined low mass bounds 
from [58], that are dominated by search channels where the Higgs is produced in associa- 
tion with a vector boson and decays into bb, are not strong enough yet to exclude BMSSM 
parameter space. 

In Fig. 6 we present the result of a parameter scan of the model as defined in the figure 
caption, for the mass of the lightest Higgs boson as a function of the mass of the charged 
Higgs on the left and as a function of the phase of a on the right. All the shown points 
satisfy the requirement of a correct EWSB (see Sec. 2.2). In addition, we impose constraints 
from ED Ms (points in blue are excluded) , then we check the compatibility with LEP using 
Higgsbounds (points in orange are excluded) and finally we impose the Tevatron bounds 
as described previously (points in red are excluded). It is interesting to note that points 
satisfying all the constraints (in green/light gray) arise in the entire range for the phase of 
a, i.e. in spite of the strong constraints coming from EDMs, large CP violating phases are 
allowed in the model. 

As expected, due to the higher dimensional operators, the values of the lightest Higgs 
mass cover a much larger range as compared to the MSSM and reach up to ~ 210 GeV, 
even in presence of large CP violating phases (see right panel of Fig. 6). As discussed in [19] 
in the CP conserving case, effects of higher dimensional operators at the 1/M^ order can 
increase that limit further up to ~ 300 GeV. As we will discuss in the next section, a 
rather heavy lightest Higgs boson (in the range (170-210) GeV) has usually a gluon gluon 
fusion production cross section times branching ratio into vector bosons that is enhanced 
with respect to the SM. Such a Higgs boson will be easily probed at the LHC already with 
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1 fb~^ (see Fig. 11 below). 

The few allowed points below the SM LEP bound of 114.4 GeV, correspond to relatively 
light charged Higgs masses (~ (130 — 170) GeV). These points can be reached both in the 
CP conserving and in the CP violating cases (see right panel of Fig. 6) and arise once 
the lightest Higgs boson is close to a pseudo scalar state, namely with a strongly reduced 
coupling to gauge bosons. Correspondingly, such Higgs bosons are very difficult to probe 
both at the Tevatron and at the LHC. 

4.2 Characteristic Scenarios for Collider Searches 

In this section we present our analysis of the Higgs collider phenomenology of the model. 
Taking into account constraints from vacuum stability and EDMs as well as collider con- 
straints from LEP and the Tevatron we concentrate on those aspects that can distinguish 
the BMSSM with CP violation from the case without CP violation that has been exten- 
sively studied in [27] . 

Genuine signatures of the new CP violating phases of the BMSSM can in principle 
occur in CP asymmetries based on the longitudinal r polarization in the WW -^ Hi ^- tt 
channel. Such asymmetries can be resonantly enhanced in scenarios where two (or all three) 
neutral Higgs bosons are nearly degenerate, with mass differences comparable to their decay 
widths [59] . A measurement of these observables however appears to be challenging at the 
LHC and a detailed study of them is beyond the scope of the present work. We concentrate 
instead on distinct features in the Higgs spectrum, the Higgs - vector boson couplings, the 
gg — )• Hi production cross sections and the Higgs branching fractions that we analyze first 
in a generic scan of the parameter space (Sec. 4.2.1) and secondly in several representative 
scenarios (Sec. 4.2.3). We calculate the Higgs branching ratios using the results collected 
in [51] in combination with HDECAY [54], after implementing the BMSSM Higgs spectrum 
and couplings appropriately. 

As the higher dimensional operators of the BMSSM have the strongest impact on the 
mass of the lightest neutral Higgs boson for small values of tan/3, we concentrate our 
discussion here to scenarios in the low tan/3 regime. This has the additional advantage 
that constraints from EDMs are kept at a minimum. 

4.2.1 Generic Features of the Parameter Scan 

In the mass range M//. > 140(120) GeV, the process gg ^ Hi ^ WW/ZZ is the dominant 
one for Higgs searches at the Tevatron (LHC). For lower Higgs masses, Mh^ < 125 GeV, 
the associated production with subsequent decay into bb, V — )• VHi — )• bb, is the main 
search channel at the Tevatron. At the LHC on the other hand, the gg ^ Hi ^ 77 
channel is the most important one for low Higgs masses. The production in vector boson 
fusion and decay into t^t^, WW -^ Hi ^ t^t^, as well as the associated production 
of boosted Higgs bosons that decay into bb, are only marginally important at the 7 TeV 
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LHC run. We remark that the BMSSM effects in the Higgs couphngs can of course modify 
the relative importance of the different search channels. In the following we consider a 
parameter scan, analyzing the most relevant Higgs search channels mentioned above. We 
use the same parameter scan already described in the previous section for the study of 
LEP and Tevatron bounds. 

The left panels of Fig. 7 show the product of the gg — )• Hi production cross sections and 
Higgs branching ratios into WW, normalized to the SM values. In the right panels instead, 
the cross section of the complementary LHC channel WW -^ Hi ^ tt normalized to the 
SM cross section is shown. For the low values of tan /? considered here, the effective Hibh 
and HiTf couplings are the same to an excellent approximation. In addition, also the ratios 
of the Higgs production cross sections in vector boson fusion and associated production 
with the corresponding SM values are identical. Correspondingly, the WW -^ Hi ^ tt 
plots are also valid for the low mass Higgs search at the Tevatron in the V — )■ VHi — )• Vbh 
channel and similarly also for the boosted Higgs analysis at LHC [60, 61]. In Fig. 7, we do 
not present instead the cross sections of the gg ^ Hi ^^ 77 channel, since we find in the 
case of CP violation, that the lightest Higgs boson generically does not show any sizable 
enhancement compared to the SM predictions (see also Sec. 4.2.2 below). The two heavier 
Higgs bosons show a considerable enhancement only for rather heavy masses (> 250 GeV) 
for which the di-photon branching ratio is tiny anyway. 

From the first row of Fig. 7, we observe that typically there are good prospects for 
the detection of the lightest Higgs boson in the WW channel, both in the CP conserving 
case (orange/dark gray) and in the CP violating case (green/light gray). For Mh^ ^ 
140 GeV, the product a {gg ^ Hi) BR(//i — )• WW) is suppressed with respect to the 
SM; still for many points of the scan, the suppression is not large (in the range 0.8- 
0.6) and does not prevent the possibility to probe the lightest Higgs boson in the WW 
channel at the Tevatron and the LHC. We also note that the upper bound on the product 
^ {99 ~^ Hi) BR [Hi — )• WW) is due to the maximal value for the charged Higgs mass that 
we allow in the scan: taking larger masses for the charged Higgs boson would allow to 
reach values closer to 1 for smaller Mh^ . 

Concerning the decay of the lightest Higgs to tt, we note typically a slight enhancement 
over the SM prediction, even in the low mass region ((115-130) GeV). This enhancement is 
connected to a suppression of the Hi — )• WW branching ratio as observed above. Given the 
rather limited sensitivity of ATLAS and CMS for the tt channel [62, 63], the enhancement 
is however not sufficient to probe the Higgs boson in this channel alone at the ongoing 
7 TeV LHC run. 

Similarly, also the two heavier Higgs bosons can only be discovered in the WW channel. 
The TT channel is in fact sensibly enhanced over the SM prediction only for masses above 
(140-150) GeV, for which the (not SM-normahzed) a {WW -^ i72,3) BR (ifa.s ^ tt) is too 
small to allow the detection of the heavier Higgs bosons. 



20 



T 

PQ 

T 



i ,\ — 


/ : 1 








m 




H^ 










1 


^^ 


^,v 










^ 


p 


f 










V" 








aR 






m 




•T 










w^ 












r-T— f 


S ;. 


^'.: 









T 

g 

ca ij 

X 

T 






Mj,, [GeV] 



SO 100 120 14[> 160 ISO 200 

M„, [QeV] 



T 

£- 0-8 

PQ 

X 0-6 



\ ( ; 1 




-J- 






m 






f 






■"v 






i. 


,^1 • 


v.. 




z-'' 




' '.' .' 


•'•> ■•..• 






iAm^ 



^ 

to 

X 0-6 

T 



M„^ [GeV] 



1^ K *- 



ca 

X 



T 

"b 






-jitS./i 



5*^ 



&i 



>■■'■ 



■ " •■:l-V:'£^^4vvp^-';i'"-:i;.;.:"-.^ 



M^SMMM 



Mh, [GeV] 




I/, |(t^\ I 



M;,3 [GeV] 



Figure 7: Left: gg — ^ i/^ production cross sections times Higgs branching ratios into WW nor- 
malized to the SM values. Right: WW — ^ Hi production cross sections and Higgs branching 
ratios into tt normalized to the SM values. Orange (dark gray) points correspond to approximate 
CP conservation (|Arg(a,Cij)| < 0.1); green (light gray) points correspond to the CP violating case 
(|Arg(a, Lo)\ > 0.1). Only points allowed by LEP and Tevatron bounds, vacuum stability and EDMs 
are shown. The gray region is excluded by the latest combined Tevatron analysis with 8.2 fb~^ [56]. 
See caption of Fig. 6 for the details of the scan. 
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In the case of CP conservation, the second Hghtest Higgs is typicahy the pseudoscalar 
state and therefore does not couple to gauge bosons. In the presence of CP violation 
however, also the second lightest Higgs has non zero couplings to WW and ZZ and sizable 
55 — >• i^2 — >• WW and WW -^ H2 ^ tt cross sections arise generically. On the other hand, 
it is interesting to note that the gg ^ H2 ^ WW cross section of the second lightest Higgs 
can even be larger than the SM prediction both in the CP violating case and even in the 
CP conserving case. This happens in regions of parameter space where the pseudoscalar is 
either the lightest or the heaviest Higgs. 

In general we observe, that from the scatter plots shown in Fig. 7 it is rather difficult 
to distinguish the CP conserving from the CP violating case, in particular for the lightest 
and the heaviest Higgs boson, since many features that are characteristic for CP violation 
are obscured by the general parameter scan. In the following we therefore discuss several 
concrete scenarios that allow a distinction between the CP violating and the CP conserving 
BMSSM. We analyze scenarios with an enhanced di-photon signal, scenarios with Higgs 
cascade decays and in particular scenarios with three neutral heavy Higgs bosons that show 
features that are characteristic for the presence of CP violation. 

4.2.2 CP Conserving BMSSM Scenarios 

The BMSSM without CP phases is of course a subset of the parameter space that we are 
considering here. In the following we comment on two benchmark scenarios which are 
present already in the absence of CP phases but become rather fine tuned once CP phases 
are turned on. 

Scenarios ^vith Enhanced Di-Photon Signal: As discussed in [27] for the framework 
of the BMSSM without additional sources of CP violation, scenarios with an enhanced 
di-photon signal are in principle possible only in regions of parameter space where the 
Hibb coupling is strongly suppressed with respect to the SM. In that way enhancements 
in the gg — )• Hi — )• 77 cross sections up to an order of magnitude compared to the SM can 
be achieved. Such scenarios are also present in our framework. However, for non-trivial 
phases they are much more difficult to realize. Neglecting for simplicity tan/3 resummation 
corrections, the Hibb coupling in presence of CP violation reads 



C = tan/3 

V 



^02i--Ou]bbHi-03ibn5bHi . (39) 

•S/3 S/3 



In the absence of CP violation one has O21 = O31 = 0, On = 1 and a suppressed Hibb 
coupling simply corresponds to Sa "^ c^. In the presence of CP violation however one has 
to simultaneously ensure O31 ^ l/tan/3 and Ca02i — SaOii ^ cp. This requires a large 
amount of fine tuning between several parameters. We conclude that a strongly enhanced 
di-photon signal would not be a characteristic signal of CP phases in the BMSSM, but 
would be a hint towards the absence of additional sources of CP violation. 
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Scenarios with Higgs Cascade Decays: As discussed in [27, 64], multi-Higgs decay 
chains can be reahzed in the BMSSM without CP phases. The scenarios considered in 
these works contain a (very) hght pseudoscalar Higgs that avoids the LEP bound as it 
does not couple to the Z boson. The higher dimensional operators in the BMSSM ensure 
that the two scalar Higgs bosons are heavy enough in order not to violate LEP constraints. 
The h/H — )• AA — )• 262r or At decays can then be the most promising search channels for 
the scalar Higgs bosons. However, in the presence of CP violation the three Higgs bosons 
mix and therefore the would-be pseudoscalar will always couple to some extent to gauge 
bosons. Therefore, for non-trivial phases it is typically more difficult to have a Higgs boson 
considerably below 114.4 GeV while avoiding the LEP bounds. 

In addition we stress that in scenarios with a light pseudoscalar-like Higgs boson, it 
is very difficult to simultaneously fulfill constraints coming from vacuum stability and 
the lightest chargino mass. As seen in Sec. 2.2, the vacuum stability constraint then 
becomes even stronger once CP phases are introduced. Furthermore, due to the very 
light Higgs spectrum also EDM constraints are extremely stringent in such a scenario and 
allow for non-trivial phases only in very fine tuned strips in parameter space. Again we 
conclude that scenarios with Higgs cascade decays into a light pseudoscalar-like Higgs are 
not characteristic of the CP violating BMSSM but more generic in absence of CP violation. 

4.2.3 Characteristic CP Violating Scenarios with 3 Heavy Higgs Bosons 

One of the distinct features of the BMSSM is the possibility to have three heavy neutral 
Higgs bosons with masses Mh^ ^ 140 GeV. In the following we discuss three concrete 
scenarios in which this is realized. The input parameters for all the scenarios are conve- 
niently collected in Tab. 1. A detailed discussion of the Tevatron and LHC sensitivities to 
the scenarios is given at the end of the section. (For a detailed analysis of the Tevatron 
sensitivity to MSSM Higgs bosons see [65, 66]) 

Scenario la: In the first scenario, all three neutral Higgs bosons are heavier than 150 GeV 
and decay dominantly into WW. This scenario cannot be realized, either in the BMSSM 
without CP violating phases or in the CP violating MSSM but is unique to the BMSSM with 
CP violation. In the BMSSM without CP violation only the two scalar Higgs bosons couple 
to the weak gauge bosons and in the MSSM with CP violation one Higgs is always light 
(Mhi ^ 130 GeV) and therefore has only a small to moderate branching ratio into WW. 
Tab. 2 shows the predictions for the Higgs spectrum, effective Higgs - ZZ and Higgs - gg 
couplings as well as the most important Higgs branching ratios for the example parameter 
set given in the first row of Tab. 1 that realizes such a scenario. Fig. 8 shows these 
quantities as functions of the phase of a, keeping the phase of oj fixed to the same ratio 
Arg(a;)/Arg(a) = —1/5 as in the first row of Tab. 1. 

In general, both in this as well as in the other scenarios considered below, one observes 
a strong dependence of the Higgs couplings and branching ratios on the phase of a, accom- 
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Sc. II 
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1 


1 


0.8 


\u>\ 


2 


1 


1.6 


Arg(a) 


7r/2 


tt/A 


-2tt/3 


Arg(a;) 


-vr/lO 


-tt/20 


tt/20 


tan/? 


2 


2 


3 


Mh± [GeV] 


195 


225 


166 


M [TeV] 


2.5 


2 


2 


fi [GeV] 


160 


190 


140 


ms [GeV] 


160 


400 


100 



Table 1: Input parameters for the scenarios discussed in the present section. For all scenarios 
we choose a common squark mass of m = 800 GeV, a common slepton mass of fhi — 1100 GeV, 
trilinear couplings At ~ 2m, Ab = Ar ~ and gaugino masses Mg = 3M2 = 6M1 = 1200 GeV. 

panied by a non-negligible dependence of the mass spectrum as well. For non trivial values 
of Arg(a) all three Higgs boson couple to the gauge bosons and all three couplings can 
be large enough that all three Higgs bosons decay dominantly into WW (see for example 
Arg(a) ~ 1.5). 

Concerning the spectrum, we remark that for — 1 ^ Arg(a) < 1, the lightest Higgs 
boson lies between 160 GeV and 165 GeV, a region where the SM Higgs is excluded by the 
Tevatron. Due to the fact that in our scenario this Higgs decays mainly into WW and has 
a production cross section that is only slightly reduced compared to the SM case, Tevatron 
data indeed excludes small values of Arg(Q) in this scenario. In addition, we observe 
that the heaviest Higgs is around 200 GeV in the full range for Arg(a) and therefore has 
a sizable branching fraction not only to WW but also to ZZ. Correspondingly, such a 
heaviest Higgs would first be detected in the Hj, — ?• ZZ — )• 4£ channel at the LHC. 

We remark however, that scenarios in which all three Higgs bosons decay dominantly 
to WW do not necessarily imply a heavy Higgs with a mass of Mn^ > 200 GeV. We also 
found scenarios with all Higgs masses in the range 155 GeV - 180 GeV, where the main 
decay channel for all Higgs bosons is Hi — )• WW. 

Scenario lb: This scenario consists of a heavy SM like Higgs with a mass of Mhi — 
150 GeV that decays mainly into WW, and two additional heavy Higgs bosons with masses 
M}{2 3 ^ 200 GeV. An example input parameter set for such a scenario can be found in 
the second row of Tab. 1. In Tab. 2 the corresponding Higgs spectrum, effective Higgs 
- ZZ and Higgs - gg couplings as well as the most important Higgs branching ratios are 
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shown. Finally, Fig. 9 shows them as functions of the phase of a, for the same ratio 
Arg(cj)/Arg(a) = —1/5 as in the second row in Tab. 1. 

From the plots in Fig. 9 we observe that the Hi — t- WW decay is indeed the main 
channel for the lightest Higgs for the full allowed region of Arg(Q;). The largest branching 
fraction for the heaviest Higgs on the other hand is always H^ — )■ bb. While the heaviest 
Higgs mainly corresponds to the heavy scalar of the theory, its coupling to gauge bosons 
is tiny. This happens because this scenario already approaches the decoupling limit where 
the lightest Higgs is SM like and carries basically the full coupling to gauge bosons. In 
absence of CP violation, the second heaviest Higgs is the pseudoscalar and therefore does 
not couple at all to gauge bosons. For non-trivial values of the phase of a however (see e.g. 
Arg(a) ~ — 1), the would-be pseudoscalar can mix sufficiently strong with the light scalar, 
acquiring a non-negligible coupling to gauge bosons and decaying mainly into WW. 

As the lightest Higgs has a mass of about (150 — 155) GeV and is SM like, it might 
be even in the reach of Tevatron. Also the second Higgs decays mainly into gauge bosons 
and the best search channel at LHC is H2 — ?• ZZ — )■ 41. Even though the heaviest Higgs 
has the largest branching ratio into bb we expect the most promising search channel to be 
H^ -^ ZZ — )• 4i, due to the non- negligible branching ratio into ZZ at the level of 10%. 
Still very high statistics will be required to probe this Higgs boson. 

Scenario II: In this scenario all three Higgs bosons have similar masses in the range 
145 GeV < MHi ^ 160 GeV and one would naively expect them to decay dominantly into 
WW. However the situation is reversed compared to Scenario la discussed above, since 
all the three Higgs bosons are decaying mainly into bb. An example parameter point that 
leads to such a scenario is given in the third row of Tab. 1. Its Higgs spectrum, effective 
Higgs - ZZ and Higgs - gg couplings as well as the most important Higgs branching ratios 
are again summarized in Tab. 2 and also shown as function of the phase of a in Fig. 10, 
with Arg(c(j) fixed to 7r/20. 

From the plots in Fig. 10 we observe that the Hi -^ bb decay is indeed the main channel 
for the three Higgs bosons in the region close to Arg(a) ~ —2. While for non trivial phases 
the HiWW coupling is again shared among all three Higgs bosons, this time it is made 
small enough such that the dominant decay mode for all three Higgs bosons is Hi — )■ bb. 
Correspondingly such a scenario is essentially unconstrained by Tevatron. 

Concerning the Higgs searches at the 7 TeV LHC, probing this scenario appears to be 
very challenging. As it is evident from Tab. 2, the gg ^ Hi ^ 77 cross sections are strongly 
suppressed compared to the SM by a factor of 5 - 15. Even though the V — )■ VHi — ;• Vtt 
channels are enhanced compared to the SM by a factor of 1.5 - 2, the large masses of the 
Higgs bosons in the considered scenario prevent the Higgs bosons to be probed in this search 
mode. A similar conclusion is expected for the boosted Higgs search. Therefore, despite 
the suppressed branching ratios, the most promising channel seems to be still Hi — )■ WW, 
even if large statistics will be required in order to be sensitive to the small signal that the 
scenario predicts. 
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Scenario la 


Hi 


H2 


Hs 


Mh, [GeV] 


157 


177 


202 


F'^ 


0.94 


0.04 


0.02 


^ggH^ 


0.72 


0.62 


0.47 


BR{Hi -^ bb) 


15% (8%) 


34% (0.6%) 


24% (0.2%) 


BR{Hi -^ WW) 


76% (83%) 


58% (95%) 


53% (74%) 


BR{Hi -^ ZZ) 


6% (7%) 


2% (4%) 


19% (26%) 


BR{Hi -^ 77) X 10^ 


9(9) 


0.8 (1.2) 


0.2 (0.5) 



Scenario lb 


Hi 


H2 


H3 


Mh, [GeV] 


153 


201 


233 


^ZZHi 


0.96 


0.03 


0.004 


F"^ 


0.84 


0.64 


0.35 


BR{Hi -^ bb) 


19% (13%) 


21% (0.2%) 


51% (0.1%) 


BR{Hi -^ WW) 


69% (74%) 


56% (74%) 


29% (71%) 


BR{Hi -^ ZZ) 


7% (8%) 


19% (26%) 


12% (29%) 


BR{Hi -^ 77) X 10^ 


12 (12) 


0.6 (0.5) 


0.5 (0.3) 



Scenario II 


Hi 


H2 


H3 


Mh, [GeV] 


147 


150 


162 


F'^ 


0.62 


0.32 


0.06 


^ggm 


0.41 


0.53 


0.39 


BR{Hi -^ bb) 


69% (22%) 


72% (16%) 


65% (2%) 


BR{Hi -^ WW) 


20% (63%) 


17% (69%) 


26% (94%) 


BR{Hi -^ ZZ) 


3% (8%) 


2% (8%) 


1% (3%) 


BR{Hi -^ -f-f) X 10^ 


6(16) 


3(13) 


0.5 (4) 



Table 2: Predictions for the Higgs spectrum, the effective ZZHi and ggHi couplings as well as 
the branching ratios for the most important decay channels of the Higgs bosons in the example 
scenarios defined in Tab. 1. The corresponding SM branching ratios for a SM Higgs with the same 
mass are given in parenthesis. In all scenarios of this section the Hi — > tt branching ratios are 
given approximately by BK{Hi -^ tt) ~ jijBR(iJi -^ bb). 
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Figure 8: Higgs spectrum, Higgs - vector boson couplings, gg — >■ Hi production cross sections and 
Higgs branching ratios as function of the phase of a for the Scenario la. The phase of w is fixed to 
Arg(a;) = — Arg(a)/5 in the plots, to keep under control EDM constraints. The side bands (light 
red shaded regions) correspond to vacua that are only local minima of the potential and the central 
band (light gray shaded region) is excluded by Tevatron Higgs searches. 
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Figure 9: Higgs spectrum, Higgs - vector boson couplings, gg — >■ Hi production cross sections and 
Higgs branching ratios as function of the phase of a for the Scenario lb. The phase of w is fixed to 
Arg(a;) = — Arg(a)/5 in the plots, to keep under control EDM constraints. The side bands (light 
red shaded regions) correspond to vacua that are only local minima of the potential. 
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Figure 10: Higgs spectrum, Higgs - vector boson couplings, gg — > Hi production cross sections 
and Higgs branching ratios as function of the phase of a for the Scenario II. The phase of w is 
fixed to Arg(a;) = 7r/20 in the plots, to keep under control EDM constraints. The inner band (light 
gray shaded region) is excluded by EDM constraints. The external bands (light red shaded regions) 
correspond to vacua that are only local minima of the potential. 
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Figure 11: Left: 95% C.L. sensitivities of Tevatron and LHC in the H -^ WW/ZZ channels. 
Right: 95% C.L. sensitivities of Tevatron in the H —>■ bb decays channels with the Higgs produced 
in association with a vector boson. The black lines in the left plot are the combined expected 
ATLAS and CMS sensitivities for 1, 5, and 10 fb^^. The gray region in both plots is excluded by 
Tevatron, the dashed gray line shows a naive extrapolation of the Tevatron sensitivity with 10 fb^^. 
The blue triangles, green diamonds and orange squares correspond to the Higgs bosons of Scenario 
la, lb and II, respectively. The red circle in the right plot corresponds to the lightest Higgs boson 
in a scenario with large Bs mixing phase to be discussed in Sec. 5. 

Fig. 11 shows the Tevatron and LHC sensitivities to the three scenarios discussed above. 
The blue triangles, green diamonds and orange squares correspond to the Higgs bosons 
of Scenario la, lb and H, respectively. The left plot shows the product of the gg — >• Hi 
production cross section times the branching ratio into vector bosons normalized to the SM 
value. The gray region is excluded at the 95% confidence level by current Tevatron data and 
is based on a naive combination of CDF and DO bounds [67, 68, 69]. We mention that the 
Hi — )• WW/ZZ results in [67, 68, 69] are not only based on gluon-gluon fusion but contain 
also vector boson fusion and associated production at a level of up to 25%. Still, as the 
gg -^ Hi ^ WW/ZZ channel dominates the obtained constraints, we consider the shown 
bound to be a reasonable estimate of the combined Tevatron exclusion reach. The gray 
dashed line shows the expected Tevatron sensitivity, extrapolating the CDF and DO bounds 
to an integrated luminosity of 10 fb^^ each. To be conservative, for the extrapolation we 
only take into account the higher statistics and do not assume any additional improvements 
of the analysis. The black lines correspond to a combination of the expected sensitivities 
from ATLAS [62] and CMS [63] at a center of mass energy of 7 TeV. The ATLAS and 
CMS curves are based on a combination of several production mechanisms. However, in 
the shown region with M//. > 140 GeV, the dominant production mechanism is gg ^ Hi, 
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with vector boson fusion and associated production contributing only ~ 10% and ~ 5%, 
respectively (see e.g. Tab. 2 in [62]). 

We observe that the lightest Higgs bosons of Scenarios la and lb are close to the current 
Tevatron exclusion and just outside the final expected Tevatron sensitivity. They should be 
easily accessible at the LHC already with an integrated luminosity of 1 fb~^. As expected, 
all three Higgs bosons of Scenario la are well within the LHC sensitivity and can be probed 
with 5 fb~^. In Scenario lb on the other hand, only the two lightest Higgs bosons can be 
probed in the Hi — )■ WW channel, while the heaviest Higgs that decays mainly into bb is 
inaccessible in the 7 TeV LHC run, unless more than 20 fb^^ per experiment are collected. 

In Scenario II finally, despite their large masses, all three Higgs bosons decay mainly 
into bb and are correspondingly harder to be probed in the Hi — t- WW channel. Combining 
the sensitivities of both ATLAS and CMS, one would need more than 10 fb~^ at the 7 TeV 
LHC to probe the three Higgs bosons in this scenario.^ 

The right plot of Fig. 11 shows the Tevatron sensitivity to the Higgs bosons of the three 
scenarios in the production in association with a vector boson and the subsequent decay 
into a bb pair. The gray region is excluded by current Tevatron data at the 95% confidence 
level and based on a naive combination of CDF and DO results [70, 71, 72, 73, 74, 75, 76, 77] . 
We checked that our combination roughly approximates the official Tevatron combination 
in the low mass region [58]. The dashed gray line shows a conservative extrapolation of 
the Tevatron sensitivity with 10 fb~^ from both CDF and DO, taking into account only the 
higher statistics in the extrapolation. Even though several Higgs bosons of our scenarios 
show a signal that is enhanced by a factor of 1.5 - 2 compared to the SM, we observe that, 
due to their rather heavy masses, they are far outside the expected Tevatron bounds. The 
same is true for the WW -^ Hi ^ t~^t~ mode and the boosted Higgs searches at the 
7 TeV LHC. 

5 The Bs Mixing Phase in the BMSSM 

In this section we study the possible impact of the new CP phases in the Higgs sector of 
the BMSSM on CP violation in meson mixing. 

The Bs mixing phase is predicted to be tiny in the Standard Model and therefore offers 
excellent opportunity to probe new sources of CP violation in NP models. Interestingly, 
combining data from CDF and DO on the time dependent CP asymmetry in Bg — )• ip(p [78, 
79, 80, 81] as well as a DO analysis of the like-sign dimuon charge asymmetry [82], one 
finds a large Bg mixing phase at the level of (2-3)(t above the SM prediction [83, 84]. 
A very recent update of the DO analysis [33] finds a discrepancy in the like-sign dimuon 



® The two lightest Higgs bosons of Scenario II are very close in mass with a mass splitting of only 
3 GeV. Given the finite mass resolution of ATLAS and CMS, these two Higgs bosons could appear as 
a single particle with the combined cross section and thus might be accessible in the particular scenario 
considered here. 
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charge asymmetry of 3.9a with respect to the SM and strengthens this resuh. In addition, 
tensions in the fits of the CKM matrix seem to hint towards sizable NP contributions to 
CP violation in Ba mixing [85, 86, 84, 87]. 

As shown in [88, 89] , generic two Higgs doublet models with Minimal Flavor Violation 
(MFV), where the CKM matrix is the only source of flavor violation [30, 31, 32] and 
additional CP violating phases in the Yukawa and the Higgs sector are allowed, are able to 
generate large effects both in the Bg and B^ mixing phases in agreement with the present 
experimental data and compatible with constraints from BR(i?s,d ~^ /^"'"A*") ^-^id EDMs. 
The same is not true in the MSSM with MFV. While 1-loop box contributions to meson 
mixing are generically small [90], sizable effects can in principle be generated at the 2-loop 
level in the large tan/3 regime by so called double Higgs penguin contributions [91, 92]. 
However, these double Higgs penguins are strongly constrained by the experimental bound 
on the branching ratio of the rare Bs — )• fJ'~^lJ-~ decay [93]. As a result one finds that 
CP violation in meson mixing remains SM like in the MSSM with MFV [94, 95]^ and 
additional sources of flavor violation are required to generate large NP effects in meson 
mixing phases.^ 

In the following we analyze to which extent this result is changed in the BMSSM. We 
consider a minimal flavor violating soft sector, i.e. no new sources of flavor violation in 
addition to the CKM matrix and investigate the impact of the higher dimensional operators 
with new sources of CP violation on B physics observables. 

5.1 Basics of Meson Mixing 

The Bq mixing amplitude is given by 

{B,\n,s\Bq) = Mf, - ^r?2 . (40) 

Assuming no NP effects in the absorptive part F^g that is dominated by SM tree level 
contributions, the effects of new short distance dynamics can be conveniently parametrized 

by 

Mf^ = (Aff2)sM + {Mf^hp = CB.e'^'" (Mf2)sM • (41) 

The Bd^s mass differences and the CP asymmetries S^Ks ^^^ S^^f, in the Bd — )• i^Ks and 
Bg — )• ip(j) decays are then given by 

AMg = 2|Mf2| = (AM,)sMCB, , (42) 

S^Ks = sin(Arg(Mi^2)) = sin(2/3 + 26 d) , (43) 



^A sizable Bs mixing phase can be possible in the so-called uplifted SUSY Higgs region [96, 97], even 
though such a framework is strongly constrained by B physics observables and [g — 2)^ [98]. 

* SUSY models that contain sources of flavor violation beyond the CKM matrix and that are capable 
of generating a sizable Ba mixing phase have been studied for example in [99, 100, 101, 95, 102, 103, 104, 
105, 106, 107]. 
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S^^ = -sin(Arg(Mf2)) = sin(2|/3,|-2e,) , (44) 

where the SM angles /3 and /3s are the phases of the CKM elements Vtd and Vts in the 
standard CKM phase convention 

Vtd = \Vtd\e-'^ , Vts = -\Vts\e-'^' ■ (45) 

While the SM predictions for the mass differences are in excellent agreement with the 
experimental determinations, the same is not true for the mixing phases. Using as input 
the CKM parameters from the NP fit in [84] one finds 

C-SM _ n. O1+0.02 r-SM _ n. n.AA+0.002 fAf.\ 

This should be compared to the measured values [108, 109] 

5;^|^= 0.67 ±0.02 , S^^ = 0.78tlil, (47) 

that differ by (2-3)o" from the SM predictions. ^ 

5.2 BMSSM Contributions to B Meson Mixing 

New Physics contributions to the Bq mixing amplitude can be encoded in the effective 
Hamiltonian 

5 3 

Hj = -^c^Oi-Y, ad, , (48) 

with the operators Oi and Oi given e.g. in [110, 111]. 

In the BMSSM with no new sources of flavor violation in addition to the CKM matrix, 
sizable contributions to the mixing amplitudes can be generated in the large tan/3 regime. 
The two most important operators that contribute to Bq mixing in such a framework read 

O2 = {bRqL){bRqL) , O4 = {bRqL){bLqR) ■ (49) 

To evaluate the corresponding contributions to the mixing amplitudes, renormalization 
group evolution [110, 111] is used to run down the corresponding Wilson coefficients C4 
and C2 from the matching scale to the scale of the B mesons where the hadronic operator 
matrix elements are given [112]. 

In the following we discuss the contributions to the Wilson coefficients, showing the 
leading terms in tan/3 and performing a v/m expansion. However, in our numerical anal- 
ysis we implement these contributions following the procedure described in [92] that goes 
beyond these approximations. 



^Note that the result for S^,j, from [109] does not include yet the updates on the time dependent CP 
asymmetries in Bs — ^ V"/" from CDF and DO [80, 81]. As these updates are in better agreement with the 
SM prediction, including them is expected to slightly decrease the central value for S^^. 
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Figure 12: Diagrams a) and b) show the most important double Higgs penguin con- 
tributions to Bg mixing in our framework. Diagram c) is the dominant contribution to 
Bg — )• /i"*"//". The leading tan/3 enhanced contributions to the mixing and decay am- 
plitudes come from flavor changing self-energy corrections in diagrams where the Higgs 
propagators are attached to the external quark legs. 

The Wilson coefficients of the above operators are mainly generated by double Higgs 
penguin contributions. In the case of Bs mixing, the dominant Higgsino loop diagrams are 
shown in Fig. 12 and can be approximated by 
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(51) 



These expressions assume a common mass fh for the squarks. MFV frameworks in principle 
allow for a mass splitting between the first two and the third generation of left-handed 
squarks that can be induced radiatively through the Yukawa couplings. In such a case, 
in addition to the Higgsino contributions shown here, also gluino contributions to FCNC 
processes arise and can become relevant [113, 114]. They are consistently included in our 
numerical analysis. The tan/3 resummation factors appearing in (50) and (51) can be 
approximated by 
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with the gluino, Higgsino and Wino contributions given already in Eq. (27). The loop 
function /i is given in appendix C. 

Expanding the Higgs propagators in (50) and (51) in M'^/M\, cot/3 and 1/M we find 
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in agreement with [115] 
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Some comments are in order: the contribution to d is proportional to mgrnh for Bg mixing 
and to mdiTib in case of Bfi mixing and therefore completely negligible in the latter case. 
The phase of the contribution to C4 is given by the same tiny phase appearing in the 
SM contribution 2Arg(VJ*) = —2/3^ ~ 2°, even for complex fj, and At. This is true in the 
approximation that enters the above equation that the squarks are all degenerate with a 
mass m}^ Compared to the MSSM, the contributions to C4 receive only corrections at the 
1/M level. 

The contribution to C2 is proportional to m^ both for Bg and Bd mixing, leading to 
NP effects in Bg and B^ mixing that are comparable in size. In the absence of the higher 
dimensional operators in the Higgs sector, i.e. in the MSSM, the contribution to C2 is 
highly suppressed by cot^ /3M^/M^ and, despite its enhancement by rrib/rnq compared to 
C4, it is usually completely irrelevant [92, 93, 115]. In the BMSSM on the other hand, once 
1/M effects are taken into account, C2 has the same tan^/3 dependence as C4 and is only 
suppressed by 1/M and v"^ /M\. Correspondingly, due to its enhancement by mh/nig, it is 
typically much more important than C4, in particular for small Higgs masses. In addition, 
the contribution to C2 is very sensitive to NP sources of CP violation. Its phase, that in 
turn generates a phase of the mixing amplitude M12 through (41), can be induced by ^J-At, 
the phase of the Higgs VEV 9, or directly through the Higgs propagator by the parameters 
A5 and Ag. A sizable complex C2 is the main qualitative difference between the BMSSM 
and the MSSM in the flavor sector. 

5.3 Constraints from other Flavor Observables 

There are several severe flavor constraints on MFV frameworks in the large tan/3 regime 
[117, 98]. In our case, the most important ones are the rare Bg — )• fJ.~^^~ and B — )■ Xgj 
decays. While also the B — t- tv, B — t- Dtu and K ^ iiv decays are known to be very 
sensitive probes of extended Higgs sectors, they are much less constraining than Bg — )• 



^"For large splittings between the strange and the bottom squark masses also C4 can in principle depend 
on NP phases [93, 116, 97]. 
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;U+/U and B — )• Xg'j in regions of parameter space where the Bg mixing amphtude receives 
sizable NP contributions.^^ 

5.3.1 Bs -^ fx+fJ-- 

Recently, the CDF collaboration presented the first two sided bound on the branching ratio 
of the rare decay Bg — )• Ai"*"/^" [34] . They find 

BR{B, ^ ^+Ai-)exp = (l-Stli) X 10-^ > (55) 

with a central value that is roughly a factor of six above the SM prediction [119] 

BR{Bs -^ ^+Ai")sM = (3.2 ± 0.2) x 10"^ . (56) 

CDF also provides an upper bound on the branching ratio at the 95% C.L. [34] 

BR(S^^^+/i-)cxp <4.0x 10"^ . (57) 

In our numerical analysis in Sec. 5.4 we will use the upper bound (57), and comment on 
the implications if the central value in (55) will be confirmed with larger significance. 

As it is well know, the Bg — )• fJ,^fJ-~ decay constitutes a very important constraint of 
the MSSM in the large tan/3 regime, as its branching ratio grows with tan^/3 [120, 121] 
and enhancements by orders of magnitude compared to the SM prediction are possible. 
Approximately one has 

tiK(Bs -)• fi^n JSM 
The dominant Higgsino contributions to S and P are illustrated in Fig. 12 and given by 

qH _ 2 *l rrit Mte*''/i(x^) 
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V j^ (^02i - "-^Ou - iOs^) (^-^02i - '-^0^^ , (59) 






Xl^T7^{03i + ^-02^-^-Ol^]03i. (60) 



^^ While there is a tension at the (2 — 3)a level between the experimental determination of BR(_B — > tu) 
and the low value of V^b that is preferred by fits of the Unitarity Triangle [118], a more conservative SM 
prediction of the BR(_B — > ti/), that is based on the PDG value of Vut, still leaves sizable room for NP 
contributions to B — > ti/. 
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For the SM loop function one has ^0(2;^) ~ 0.96. Gluino contributions to S and P are 
possible within MFV for a splitting between the first two and the third generation of left- 
handed squarks [114] and are included in our numerical analysis. In the decoupling limit 
and at leading order in the 1/M expansion, the Higgs propagators entering in the above 
expressions reduce to 

and one recovers to first order the well known expressions in the MSSM with MFV. The 
expressions for 5* and P in (59) and (60) scale as tan^ I3/M\ at the amplitude level. As 
they depend on the same combination of SUSY parameters as the double Higgs penguin 
contributions to Bg mixing, the NP contributions to Bs — )■ /ti "*"/"" and Bs mixing are 
strongly correlated in the large tan/3 regime. 

We stress that relative to the MSSM expectations for the branching ratio of the 
Bs — )• /^"'"Ai" decay, the BMSSM physics only leads to 1/M suppressed corrections, that 
we however include in our numerical analysis. Analogously, also the B^ — t- iJi^ iJi~ decay 
remains to first approximation MSSM-like in the BMSSM and, given the current experi- 
mental bound on its branching ratio BR(i?rf — )• /.i^/i^) < 6.0 x 10^^ [34], it is much less 
constraining than the Bg — t- h^ fi~ decay in models with MFV. 

5.3.2 B -^ Xsj 

Combining the latest experimental results for the branching ratio of the B — )■ Xgj decay 
from [108] with the NNLO SM prediction [122] (see also [123] and [124]) one finds 

BR(B -^ Xg-y) , , 

^^^7 = T^T^/R^y ^ = 1-13 ± 0.12 . 63 

In frameworks with Minimal Flavor Violation, the prediction for the branching ratio can 
be approximated by [125, 126] 

-Kfes7 — i + a77|07 I +a88|(-'8 ] 

+Re {ar Cf'') + Re (ag Cf^) + Re {arsC^^C; ^P) , (64) 

where the coefficients fij are given by 07 = —2.41 + 0.21i, as = —0.75 — 0.19i, 077 = 1.59, 
ass = 0.26 and 078 = 0.82-0.30i [126]. The NP contributions to the Wilson coefficients Cj 
and Cs are evaluated at a scale of 160 GeV in the above equation. The Wilson coefficients 
entering (64) receive the most important contributions from charged Higgs - top loops 

37 



and Higgsino - stop loops in the scenarios that we consider. Apart from small cot/3/M 
suppressed corrections that are included in our numerical analysis, they have the same 
form as in the MSSM. At the matching scale they read 

^± = l-eStan^ 

^'^ l + e^tan/3 '^^^ ' ^ ' 



jH ^ mt tan/3 ^t/^ jg 
''7,8 ^2 I -j-g* tan/3 m^ 



'-'7,8 = :=:2^r7~rT-— 53:2"^ /7,8(a;^) , (66) 



with yt = m^/MJj^ and the loop functions hj^s and /y^g are collected in appendix C. Gluino 
contributions to Cj^s in the context of MFV are discussed in [127, 114] and included in our 
numerical analysis. As long as possible phases in the tan (3 resummation factors do not play 
an important role, the charged Higgs contributions interfere constructively with the SM 
contribution. They can be partially canceled by the tan /? enhanced chargino contributions 
if sign{fiAt) = + 1 . We note the appearance of the phase of the Higgs VEV in the chargino 
contribution that, as already mentioned in Sec. 3, is connected to the tan/3 enhancement. 

5.4 Numerical Analysis 

Due to the strong correlation between the Higgs penguin contributions to Bg — )• l^^ fJ-~ 
and the double Higgs penguin contributions to the Bg mixing amplitude, one expects that 
possible NP effects in Bg mixing are severely constrained by the data on BR(i?s — )■ (jl'^ (j,~). 
We stress however the different behavior of the contributions to the Bg mixing amplitude 
and Bg — )• //"^/x^ with the Higgs mass and tan/3. In the decoupling limit one has 

I ~ v"^ 1 

C4 ex tan B — ^ , C2 oc tan B — -r and 5, P oc tan B — tt . (67) 

M\ M\ M\ ^ ' 

Correspondingly, to keep the BR(i?s — )• fJ-^fJ') under control and simultaneously keep a 
sizable C2, low values of the Higgs masses and moderate values of tan/3 appear to be the 
most promising region of parameter space. 

Making tan/3 as low as possible while keeping a sizable Bg mixing amplitude requires 
certain choices for the remaining parameters entering the expression of C2. Apart from 
the obvious requirements of small Higgs masses, a large A5 and At, also a sizable negative 
H term increases the size of C2 significantly through the tan P resummation factors. 

The requirement of a large fi term implies that in the regions of parameter space where 
a sizable Bg mixing phase is possible, the vacuum is generically not absolutely stable (see 
discussion in Sec. 2.2 and in particular Fig. 1). In the minimal BMSSM setup that we 
considered up to now, typically a very deep second minimum of the Higgs potential arises 
along the D-flat direction corresponding to a VEV of f ~ ^J ixMjuj. Using simple analytic 
expressions for the bounce action [36] to estimate the tunneling rate from the electro- weak 
vacuum to the deeper vacuum, we find that the life times of these electro- weak vacuae are 
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typically much shorter than the age of the universe. This strongly constrains the allowed 
parameter space and prevents to large extent a large Bg mixing phase. 

However we remark the following: As discussed also in Sec. 2.2 this second minimum 
arises due to negative quartic Higgs couplings that are stabilized by the 1/M^ suppressed 
dimension 6 terms (5) that in turn are induced by the dimension 5 operator in the super- 
potential. At the 1/M^ level there are however several other operators that can modify 
the Higgs potential [19, 20] and that are not even necessarily suppressed by the same scale 
M. These operators can induce additional dimension 6 terms in the Higgs potential that 
have strong impact on the vacuum structure for large field values. Adding to the Higgs 
potential dimension 6 terms as in (5) with order one coefficients, we checked that they 
can indeed remove the second minimum as long as their scale is not too large, around 
M ~ 2 TeV. Generically, the additional 1/M^ operators also induce corrections to the 
quartic couplings of the Higgs potential and therefore can affect also the Higgs spectrum 
in the usual electro-weak vacuum with v = 246 GeV. However, these corrections of the 
quartic couplings depend on combinations of coefficients of the 1/M^ operators that are 
independent of the coefficients entering the induced dimension 6 terms. Correspondingly, 
in the adopted effective theory approach, one is free to choose the 1/M^ corrections to the 
Higgs potential such that their impact on the quartic couplings is small and simultaneously 
the electro-weak vacuum is stabilized. 

The very light Higgs spectrum together with the large negative /i term also leads to 
very large charged Higgs contributions to the 6 — )■ 57 amplitude. Still, one can be in agree- 
ment with the experimental constraint due to the chargino contributions that can (at least 
partially) cancel the Higgs contribution on condition that sign{fj,At) = +1. This condition 
has the advantage that gluino and chargino contributions to the tan /3 resummation factors 
add up constructively. Additionally, the mass of the stops should be rather low. We chose 
to split the masses of the third generation of squarks from the first 2 generations. This 
opens up the possibility to have also gluino contributions to the FCNC amplitudes which 
can further enhance C2 slightly. 

In Fig. 13 we show an example scenario capable of producing a large Bs mixing phase 
and being generically in agreement with the experimental constraints on Mhi, AM^, 
BR{Bs -^ /u+^-) and BR{B -^ Xs-f). We fix |a;| = 0.4, \au}\ = 2, Arg(w) = -0.75, 
Arg(Q;) = —2, fi = —950 GeV, ms = 1000 GeV, M = 6 TeV, 3rd generation squark soft 
masses m^ = m^ = 500 GeV, 1st and 2nd generation squark masses rh = 4 TeV, slepton 
masses rhi = 4 TeV, trilinear couplings At = —2.5m^, gaugino masses Mi = 200 GeV, 
M2 = 400 GeV, Ms = 1200 GeV. We plot contours of constant S^^ in the Mh± - tan/3 
plane. 

As expected, the largest values for S"^,^ that are in agreement with the BR(i?s — )• /i+/i^) 
bound are obtained for a rather light Higgs spectrum {Mjj± ~ 240 GeV) and moderate 
values of tan/3 ~ 11. Values for the Bg mixing phase up to S^^j, ~ 0.3 can be reached in 
this example scenario. 

The chosen large phases for w and a together with the light Higgs spectrum and the 
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Figure 13: Possible values of S^^ in the Mjj± - tan/3 plane for the example scenario 
described in the text. The dark gray region is excluded by the constraints from AMs 
(dotted black line), from BR(i3 — )■ Xgj) (dashed-dotted black line) and the LEP bound on 
the lightest Higgs mass (dashed black line). The light gray region is excluded by the data 
on BR(i?s — )• ;U+/i^) (solid black line) only. 

not very small tan/3 will generically also induce huge contributions to EDMs. In Fig. 14 
we show the value of 5^^ in the Arg(a) - Arg(a;) plane, fixing tan/3 and Mjj± to 11 and 
240 GeV, respectively and all the other parameters chosen as in the scenario presented 
above. While in the left plot only the constraints from AMg, the lightest Higgs mass, 
BR{B — )• Xg^) and BR(Ss — )• /u"''/i~) are applied, the right plot shows the regions that 
are excluded by the EDMs. We observe that even though the EDMs severely restrict the 
allowed values for Arg(a) and Arg(a;), there are regions in parameter space with sizable S^^ 
that pass the constraints from the EDMs. This happens because the contributions to the 
Bs mixing amplitude and the contributions to the EDMs depend on different combinations 
of the phases Arg(a) and Arg{u;). 

We remark that the Thallium, Mercury and neutron EDMs in principle also depend on 
different combinations of the NP phases because of their different sensitivity to the 1- and 
2-loop contributions. For not very small tan /3 the different EDMs tend to generically give 
complementary constraints in the Arg(a) - Arg(a;) plane and largely exclude non-trivial 
phases. This does not happen in the example scenario above where we chose heavy 1st 
and 2nd generation squark as well as heavy sleptons. In that way the 1-loop contributions 
to the EDMs decouple to a large extent and the Thallium, Mercury and neutron EDM 
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Figure 14: Left: Possible values of S'^^ in the Arg(a) - Arg(a;) plane for the example 
scenario of Fig. 13 (see text for details). The dark gray region is excluded by the constraint 
from AMs, the LEP bound on the lightest Higgs mass and the experimental data on 
BR(i?5 —7- fi^n") and BR(i3 — )• Xsj). Right: Same plot as on the left, but with constraints 
from EDMs. The solid, dashed and dotted black lines correspond to the Thallium, Mercury 
and neutron EDMs respectively. 

are all dominated by very similar 2-loop contributions, leading to approximately aligned 
constraints in the Arg(a) - Arg(cj) plane. 

We also remark that the right plot of Fig. 14 reflects the situation where Arg(a) and 
Arg(a;) are the only CP violating phases of the model. If also the ^ term or the soft SUSY 
breaking parameter introduce CP violation, additional 1-loop contributions to EDMs arise. 
Large cancellations among the several contributions can then in principle occur and the 
parameter space for Arg(a) and Arg(a;) opens up. An extensive analysis of the model, also 
allowing for a complex /i term and a complex soft SUSY breaking sector is however beyond 
the scope of this work. 

Having established that there are indeed regions in parameter space where a sizable 
S^(f) is possible, we show in Fig. 15 the results of a parameter scan of the model. We keep 
the gaugino, slepton and 1st two generation squark masses as well as the scale M to the 
values of the example scenario defined above and allow the remaining parameters to vary 
in the ranges |w| < 1, \au}\ < 2, -vr < Arg{uj), Arg(a) < vr, -1000 GeV < fi < -500 GeV, 



500 GeV < 



ms 



< 1000 GeV, 500 GeV < 



mr 



mr < 1000 GeV, At 



-(2 - 3)mi 



tan/3 < 20 and Mfj± < 500 GeV. The left plot of Fig. 15 shows the correlation between 
S^^ and BK{Bs — ;■ fi^fj.^). We observe that in our framework the current experimental 
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Figure 15: Left: Correlation between 5^^ and BR(i?s — t- /U"*"//"). The solid green line 
represents the current central value of the CDF result [34] on BR(i?s — )■ /i"*"//"), while the 
dashed lines correspond to the experimental la range. The dotted line finally shows the 
95% C.L. bound. Right: Correlation between 5^,^ and S^Ks- The green areas show the fit 
results for S^^ from [84] combined with the experimental measurement of S^Ks ^^ the 1, 2 
and 2>a level. Red (light gray) points are allowed by all constraints, (dark) gray points are 
excluded by BR(i?s — t- fi^fj,^). In both plots the black point shows the central SM values. 
In the right plot the vertical black lines indicates the la uncertainty of the SM prediction 

of SjpKs ■ 



data on BR(i?s —5- fJ-^fJ-^) gives an absolute limit on the Bg mixing phase of 5^^ < 0.4^^, 
with the central value for BR(i?s — )■ n'^fJ.") leading to S^^f, < 0.25. Furthermore, for any 
given value of S^cf,, the model predicts a lower bound on BR(i?s — )• fi^fj.^). For S^^f, ~ 0.2 
for example we obtain BR(i?s — ^ /i^/-i^) > 1.5 • 10^^. Such values will be probed in the 
near future by LHCb [128]. 

The right plot shows the strong correlation between 5^,^ and S^Ks ^^ the studied model. 
Large positive effects in S^^ always imply sizable negative effects in S^Ks ^^^ '^i^e versa, 
as it is expected because of their origin from C2. The NP effects in the Bg mixing phase 
soften the tension between the SM prediction and the value for 5"^^ preferred by recent 
fits [83, 84] and simultaneously the effects in S^Ks can lead to a very good agreement with 
its measurement. 



We end this section by listing further predictions in the region of parameter space with 



^^However, if we assume the absence of the additional operators arising at the 1/M^ level that can 
stabilize the EW minimum, then the requirement of an absolute stable EW minimum implies S^^ ^0.1. 
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largest Bg mixing phase. As it is evident from Fig. 13, the mass of the hghtest Higgs is 
close to the LEP bound of 114.4 GeV while the two heavier Higgs bosons have masses 
of Mff^ 3 — (200 — 300) GeV. Despite the rather large mass splitting, the lightest Higgs 
does contain non-negligible components of the non-SM like Higgs bosons in that region of 
parameter space and therefore has enhanced couplings to tt and in particular to bb, due 
to the tan/3 resummation factors. Correspondingly, its branching ratio to 77 is strongly 
suppressed, such that observing this light Higgs boson at the LHC will be very challenging. 
On the other hand, due to the enhanced branching ratio into bb it might be possible to 
probe this Higgs in the HV — )• bb channel, analyzing the full Tevatron data, as shown in 
Fig. 11. Also the two heavy Higgs bosons have enhanced couplings to bottoms and taus 
and they can be probed at the LHC in the inclusive -ff2.3 — ^ tt channels. The stops in 
this scenario are rather light with masses m^ < 500 GeV while the squarks of the first 2 
generations and also the sleptons are expected to be very heavy in the multi TeV range. 

All EDMs are generically predicted to be close to the current experimental bounds. 
However, no absolute lower bound can be put on the EDMs due to the possibility of large 
cancellations among different contributions. 

Finally we mention that in SUSY frameworks where the CKM matrix is (effectively) 
the only source of flavor violation but additional CP phases are considered, visible effects in 
observables that are sensitive to CP violation in the 6 — )• S7 transition are also generically 
predicted [94, 50]. Examples of such observables include the CP asymmetry in the B — )• Xsj 
decay [129, 130, 131], time dependent CP asymmetries in the B — t- 4>Ks and B — )■ r]'Ks 
decays [132] as well as CP asymmetries in the B — )• K*i~^i~ decay [133, 134]. However, 
in general no clear correlation with the Bg mixing phase can be established, as these 
observables have different dependence on the SUSY parameters. 

6 Conclusions 

In addition to the MSSM particle content, new degrees of freedom beyond the MSSM 
(BMSSM) might be present slightly above the TeV scale. In an effective field theory 
approach, the leading corrections from the BMSSM physics to the MSSM Higgs sector can 
be described by two dimension 5 operators. These two operators are a potential source of 
CP violation at the tree level. In this work we studied the impact of these new sources 
of CP violation beyond the MSSM on Higgs searches and B meson observables taking 
into account constraints from Electric Dipole Moments. In the first part of the work we 
concentrated on the Higgs collider phenomenology that is specific of the BMSSM with CP 
violation. In the second part we analyzed a complementary region of parameter space 
where interesting effects in flavor physics are possible. 

In contrast to the MSSM, the Higgs sector of the BMSSM generically is CP violating 
already at the tree level. Correspondingly the physical spectrum consists of one charged 
Higgs and three neutral Higgs bosons that are mixtures of the neutral scalar and pseu- 
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doscalar states. Taking into account constraints from vacuum stability and Electric Dipole 
Moments as well as collider constraints from LEP and Tevatron, we worked out the Higgs 
collider phenomenology of the model. We discussed the Higgs spectrum and couplings as 
well as the dominant production cross sections and the Higgs branching fractions in several 
scenarios that show distinct features of the BMSSM with CP violation. 

As long as tan/3 (as well as the BMSSM scale M) is not too large (tan/3 < 5), the 
higher dimensional operators can significantly enhance the mass of the lightest neutral 
Higgs boson above the bound of Mh < 135 GeV that holds in the MSSM. As also EDM 
constraints become extremely strong for large values of tan/3, we concentrated on the low 
tan /3 regime in our analysis of the Higgs collider phenomenology. The CP violating phases 
of the higher dimensional operators allow for sizable couplings of all three Higgs bosons to 
the weak gauge bosons. Correspondingly, we found that the most striking scenario that 
can be realized in the BMSSM with CP violation are three neutral Higgs bosons that all 
decay to the WW final state. Such a scenario cannot be realized, either in the MSSM 
with CP violation or in the BMSSM without CP violation and is unique to the framework 
considered in this work. It can be probed at the 7 TeV run of LHC. 

We also find that the model allows for benchmark scenarios where all three Higgs 
bosons are heavy with masses M//- > 150 GeV but still all decay dominantly into bb. 
Correspondingly, it will be challenging to probe such scenarios at the 7 TeV run of LHC 
as they predict signals in the studied Higgs search channels that are at the very border or 
even below the expected sensitivities with 10 fb~^. 

We also discussed distinct signals of the modified Higgs sector on flavor observables in a 
complementary region of parameter space, where the direct Higgs searches do not signifi- 
cantly depart from the MSSM expectations. In the large tan /3 regime, significant deviations 
from the MSSM predictions can in principle be expected in the Bd^s -^ 1^^ 1^^ decays and 
in Bd^s mixing. Such processes can be affected by NP contributions from Higgs and double 
Higgs penguins and are therefore highly sensitive to the Higgs spectrum and couplings. We 
find that the main qualitative difference with respect to the MSSM are NP contributions 
to the AF = 2 operator (bPLqf' . In the MSSM this operator is highly suppressed by 
cot^ /3 My^/M\ and completely negligible. Instead, the leading double Higgs contributions 
to the meson mixing amplitudes in the MSSM generate the operator {hPiiq)(bPLq). As- 
suming Minimal Flavor Violation (MFV), i.e. the absence of new sources of flavor violation 
in addition to the CKM matrix, these contributions are proportional to mhiris and mi,md 
in the case of Bg and Bd mixing respectively. Such contributions are strongly constrained 
by the upper bound on the branching ratio of the Bg — )• A*"*"/^" decay in the MSSM with 
MFV and therefore CP violation in meson mixing remains SM like. Correspondingly the 
hints from Tevatron towards a large Bs mixing phase cannot be addressed in the MSSM 
with MFV. 

In the BMSSM with MFV on the other hand, the (bPLq)'^ operator typically gives 
the dominant contribution. It is proportional to m^ both for Bg and B^ mixing and 
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can induce large CP phases to both mixing ampHtudes that are equal in size and phase. 
Even though the Bs — )■ ^^/U~ decay still gives the strongest constraint on double Higgs 
penguin contributions to meson mixing, we find that a Bs mixing phase up to S^^ < 0.4 
can be achieved in the BMSSM with MFV, compatible with EDM constraints. We find 
that stability of the electro-weak vacuum in the region of parameter space with a large 
Bg mixing phase requires the presence of additional 1/M^ operators with appropriately 
chosen coefficients. In addition, a large Bs mixing phase of S^^ ~ 0.4 implies a sizable 
suppression of the SM prediction for the B^ mixing phase S^Ks ' which is welcome in view 
of the observed tensions in the determination of the unitarity triangle. Simultaneously, also 
the BR(Ss — ^ Ai^/" ) is predicted close to the upper bound of the region recently reported 
by CDF. Interestingly, a sizable S*^^ > 0.25 implies a lower bound on BR(i?s — ^ /^^^ ) ^ 
2 X 10"*^ which is just above the central value given by CDF. 

To conclude, in this paper we have shown that the BMSSM with CP violation can 
lead to complementary benchmark scenarios that either show novel signatures in Higgs 
collider physics or can ameliorate some tensions in present B physics data. The predicted 
signals are characteristic of the presence of NP phases in the Higgs potential and allow 
to distinguish the model both from the CP violating MSSM and the BMSSM without CP 
violation and will be tested in the near future at the LHC. 

Note Added 

After the completion of this work, new Tevatron and LHC bounds on Higgs searches [135, 
136, 137] and on the branching ratio of the decay Bg — t- fi^fJ-^ [138, 139, 140] appeared. 

The ATLAS collaboration presented the combined exclusion limit for a SM Higgs with 
an integrated luminosity ranging from 1.04 and 1.21 fb~ [136]. The CMS collaboration 
presented the combination of six SM Higgs boson searches corresponding to 1.0-1.1 fb~ of 
integrated luminosity [137]. Our Higgs scenarios lie beyond the range presently probed by 
the LHC, even if the lightest Higgs boson of the first two scenarios will be soon tested by 
the LHC data. We observe that, for Higgs masses below ~ 180 GeV, the observed limits 
of ATLAS and CMS are weaker than the projected limits we presented in Fig. 11. 

Concerning the branching ratio of the rare decay Bs — )■ /i^/U~, a preliminar combination 
has been performed using the LHCb bound obtained with 337 pb^ of integrated luminosity 
(1.5 X 10~®) and the CMS bound obtained with 1.14 fb~ of integrated luminosity (1.9 x 
10"^). The combined value is BR{Bs -^ /i+zx") < 1.1 x 10"*^ at the 95% C.L. [140]. This 
result does not confirm the excess observed by CDF [34]. Using this updated bound on 
the decay of the Bg meson, we find that the possible range for the Bs — Bg mixing phase 
is S^^ < 0.15 (as observed from Fig. 15). This value is however considerably larger than 
the one possible in the MSSM with MFV. 
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A Chargino, Neutralino and Squark Masses 

In this appendix we summarize the impact of the 1/M suppressed operators in the Higgs 
sector on the chargino, neutralino and squark masses. 

We decide to factor out the phases 9u and 64 of the fields Hu and H^ (see (6)) also on 
the Higgsino fields Hu^^ i.e. we factor out the phases on the entire Higgs superfields Hu,d- 
In this way the terms that mix the Higgsinos and gauginos after electro-weak symmetry 
breaking are real and the physical Higgs phase 9 appears only in the Higgsino part of the 
chargino and neutralino mass matrices. In particular, the 2x2 chargino mass matrix reads 



^x± 



M2 



\^VCp fie'^ - UJi'^S(3Ci3 



2 



2ie 



(68) 



and the 4x4 neutralino mass matrix is given by 
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Given the non-zero phase of the Higgs VEVs we also perform a phase shift on the right- 
handed quark and lepton fields in order to have real Yukawas at the tree level. As in case of 
the Higgs fields, we apply these phase shifts on the entire quark and lepton superfields. In 
that way the left-right mixing terms in the squark and slepton mass matrices only contain 
the physical phase 9 and not 9^ and 9^ separately. The up and down squark masses read 



Mf. 



fh? + mf, 
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(71) 



where, for simplicity, we assumed a common mass fh for the left-handed squark doublets 
and the right-handed squark singlets. 
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B Electroweak Precision Constraints 



In this appendix we report the expressions for the contributions to the S and T parame- 
ter [141] from the Higgs sector of the BMSSM with CP violation. Generahzing the results 
given in [142, 19] to our case with CP violation we find 

7rM|5 = Y.^zzH,l'lk[B22{MlMJj^,MfjJ-B22iMlMl^,Ml^) 

ijk 



+ E ^''zzH, [l322{Ml Ml Ml) - B22{Ml Ml M- 

i 

- MlBoiMl Ml Ml) + MlBoiMl Ml Ml,) 



(72) 



167rs^M2,r = ^Y.^lzH\G{Ml,Ml-G{Ml,,Ml-G{Ml,M^) + G{Ml,,Ml,) 

i 

+ ^\CwHHfG{Ml,MJj^)-^ezzH,Uj,G{Ml,Ml) . (73) 



ijk 



For the loop function G one has 



G{x,y) = -{x + y) 

2 X — y 



xy IX 
log _ 

y 



(74) 



and B22 and Bq can be found in [142]. The mass Mj-ef represents a reference value for the 
Higgs boson mass in the SM. The effective couplings of the Higgs bosons with two gauge 
bosons, izzHi-, and the effective couplings of the Higgs bosons with the charged Higgs 
boson and a W, ^wHHi, that enter the above expressions, are given by 



CwHHi = Si3-a02i — CjS-aOu — iOsi , 



(75) 
(76) 



We explicitly checked that all the scenarios considered in this work are compatible with 
the constraints from S and T. 

C Loop Functions 
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